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Using series expansion by continuous unitary transformations we study the magnetic properties of 
a frustrated tetrahedral spin— | chain. Starting from the limit of isolated tetrahedra we analyze the 
evolution of the ground state energy and the elementary triplet dispersion as a function of the inter- 
tetrahedral coupling. The quantum phase diagram is evaluated and is shown to incorporate a singlet 
product, a dimer, and a Haldane phase. Comparison of our results with those from several other 
techniques, such as density matrix renormalization group, exact diagonalization and bond-operator 
theory are provided and convincing agreement is found. 
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INTRODUCTION 

Geometric frustration is a central issue in present days 
research on quantum magnets. In particular, materi- 
als with spin-tetrahedra as building blocks are under 
scrutiny, since they bear the potential of exotic mag- 
netic properties . Many investigations have focused on 
dense tetrahedral magnets, eg. on pyrochlore lattices 0. 
Quantum magnets with spatially more separated spin- 
tetrahedra have been considered far less exhaustive. The 
isostructural tellurate compounds Cu2Te20sX 2 (X=Br 
or CI) 00 are a prominent example of such 'diluted' sys- 
tems. They contain disjoint tetrahedral clusters of Cu 2+ 
with S = 1/2 0. The tetrahedra display two longer and 
four shorter edges. They form chains along the c direction 
and are separated by lone-pair cations in the a, b plane • 
Analysis of the susceptibility 0, is consistent with 
a strongly frustrated antiferromagnetic mira-tetrahedral 
exchange of Ji = 47.5(40.7)/\ along the shorter edges 
and a ratio of J2/J1 = 0.7(1) between the long and the 
short edges for X=Br(Cl). Both systems show incom- 
mensurate ordering with strongly reduced moments 
below a transition temperature of T c = 11.4(18.2)X for 
X=Br(Cl). This T c is consistent with bulk thermody- 
namics 0j H S ■ Observation of longitudinal magnons 
suggests that the Br-system could be close to quantum 
criticality 00. 

Despite the low-temperature ordering of the tellurates, 
the effects of m£er-tetrahedral coupling at various dimen- 
sionality and on an energy scale of the intra-tetrahcdral 
exchange are open issues in such magnets. In the tel- 
lurates dispersive features have been observed in inelas- 
tic neutron scattering on polycrystals 0. Two-magnon 
Raman scattering provides evidence for inter-tetrahcdral 
dispersion along the c-axis only . Considerations of the 
electronic structure however, suggested that a, 6-planar 
exchange via O-Te-0 bonds is important also |10L . 

Theoretically, inter-tetrahedral coupling in diluted 
tetrahedral magnets has been analyzed in terms of one- 
dimensional (ID) tetrahedral spin chains using bond- 
operator methods, exact diagonalization ^3] and effec- 



tive Hamiltonians [13j . The chains have been found to al- 
low for a singlet, a dimer, and a Haldane phase with low- 
energy singlet excitations close to the dimcr-to-singlct 
transition which is of first order. In 3D molecular field 
theories 00 have been used to model the magnetic tran- 
sition in the tellurates. Moreover series expansion (SE) 
up to 4th order in various 3D inter-tetrahedral couplings 
has been carried out [l4l | resulting in a rich quantum 
phase diagram with several competing states. Motivated 
by this situation, the aim of the present work is to extent 
and improve the analysis of ID tetrahedral chains carried 
out in ref. ^3] by employing the SE method of ref . ^ij ; 
however in ID only, which will enable us obtain results 
to higher orders. 

FigEl depicts the tetrahedral spin— 1/2 chain intro- 
duced in refs. [I2I fl3| . The tetrahedra correspond to 
the sites 1-4 with couplings J\ and J2 = aJ\. Inter- 
tetrahedral exchange occurs through J3 = bJ\, The 
Hilbert space of a single tetrahedron consists of 16 states, 
i.e., two singlets £1,2, three triplets 71,2,3 and one quintet 
Q whose energies are listed in Tabic [I] This level scheme 
implies that a singlet resides within the singlet-triplet gap 
of the tetrahedron for 1/2 < a < 2. The Hamiltonian of 
the chain can be written in terms of the total edge-spin 
operators Pi(2),z = s i(4) ,z + s 3(2),Zj where s^.z denotes a 
spin— 1/2 at vertex i on the tetrahedron at site I. For the 
Pi quantum numbers refer to Table [I] 
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Here L is the number of rungs (pairs of sites coupled by 
J-2 = aJi). The Hamiltonian commutes with the edge 



spin, i.e. [H, Pl =1 



(=1,2),U 



0; V/, 



1,2. Hence, the 



Hilbert space decomposes into sectors of fixed distribu- 
tions of locally conserved edge-spin eigenvalues P^;, each 
corresponding to a sequences of spin— 1 chain-segments 
{Pi, 1 = 1) intermitted by chain-segments of localized sin- 
glets (Pi t i =0). It has been shown in ref. [J^ that the 
ground state of Eq.Q occurs only within the homoge- 
neous sector of purely 52 states or in the Pi.i = 1, \/i,l 



2 




3 2 

/ / + 1 

FIG. 1: Tetrahedral cluster-chain. The unit cell, labelled 
with I, contains spin— 1/2 moments sn at vertices i = 1,.. .,4 
(solid circles). All couplings are in units of Ji (dashed lines). 
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TABLE I: Eigenstates of a single tetrahedron: singlets 
(iSi,<S2); triplets (71,72,3) and Quintet (Q). Each state is 
labelled by the Pi, 2 edge-spin quantum numbers and the en- 
ergy E/Ji. Site index / suppressed. 



SERIES EXPANSION BY CONTINUOUS 
UNITARY TRANSFORMATION 



sector. I.e., inhomogeneous phases consisting of both 
Pi 1 = and Pi 1 = 1 sites will not be ground states. Keep- 
ing 6 fixed and for a — > 00, the ground state must result 
from local S2 type of states. This is a trivial, decoupled 
singlet product phase with singlets on each rung. The 
ground state energy per tetrahedron in this phase is inde- 
pendent of 6 and it is given by: e g s 2 = — 3a/2. In the op- 
posite limit a — > the ground state is in the Py = 1, V£, I 
sector. There is a line of a first order quantum phase 
transition in the (a, 6)-plane that switches between both 
ground states, beginning at the non-interacting point 
a = 1,6 = 0. In the next Sections we will determine 
this line by SE. 

In the Pi 1 = 1, Vi, I sector the tetrahedral chain is 
equivalent to that of a dimerized spin— 1 chain with L 
sites 



-1S2 



6S 2 /S2;+i] + -L, 



(2) 



where the edge-spin operators Pi j = S2z_i and P2J = 
82; have been relabelled in terms of spin— 1 operators S;. 
For 6 = 0, the chain decomposes in isolated dimers and 
for 6 = 1 it is the homogeneous spin— 1 chain. In both 
cases the system has a spin gap. Trivially in the for- 
mer one and according to Haldane [lsl Hfil ] in the latter. 
However, for a particular b c the gap closes, defining a sec- 
ond order quantum phase transition between the dimer 
phase (6 < 6 C ) and the Haldane phase (6 > 6 C ). This 
transition will be analyzed in the next Sections by cal- 
culating the evolution of the gap in terms of the pertur- 
bative parameter 6 using SE around the non-interacting 
point (6 = 0). Our results will be compared with those 
obtained by other methods. 



The Hamiltonian of the tetrahedral chain (Eq.JI}) can 
be written as 



H 

j7 



N 



= H + aT a o + b ^2 T bi 



(3) 



-N 



Hq is the sum over local tetrahedral Hamiltonians at 
a = 0. Their spectra consist of four equidistant en- 
ergy levels E/Ji = -2,-1,0,1 (see Table HJ). With 
these levels we associate a number qi of local energy 
quanta qi = 0, ...,3. Together with Pi, 2, this charac- 
terizes the basis. Hq has an equidistant ladder spec- 
trum labelled by Q = 53; Q — refers to the un- 
perturbed ground state of Hq: |0) = . . . <Si.l/2), i-c. 
an 1S1 singlet product. The Q = 1 sector of Hq con- 
sists of linear combinations of local Ti t j triplet excita- 
tions \tj) = |iSi,i . . . 7i t j . . . with a Tij triplet on 
the tetrahedron at site j. T„ refers to the sum over the 
local terms proportional to a in Eq.Q. By construction 
this term is diagonal in the basis of Hq. The third term 
in Eq.Q refers to the inter-tetrahedral coupling via 6. 
The operators Tf, n non- locally create (destroy) n > (<) 
quanta within the ladder spectrum of Hq. For our model 
N < N max = 6, in principle. Explicit calculation of the 
Tb n however shows that N < 4 |lj| . 

Note that by a shift of one half of the unit cell, 
i.e. P2i(ii+i) Pii(2i) Eq.([Q) is symmetric under the 
operation (Ji,a, 6) — > (J16, a/6, 1/6). Therefore, in order 
to cover the complete parameter space for a, 6 > it is 
sufficient to consider the phase diagram in the range of 
a e [0, 00] and 6 G [0, 1 1. 

It has been shown [TsL flfll l2fl | that models of type 
Eq. @ allow for high-order SE using a continuous unitary 
transformation generated by the flow equation method 
of Wegner [21]]. Adapted to our case, the mapping of 
Eq.© onto the unitarily rotated effective Hamiltonian 
P e ff reads [HE^ 

Peff = Po + a T a 
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J2 b n J2 C ( m ) T bmiT b , 



1=1 |m| = n 
Af(m)=0 



(4) 



where m = (m 1 ,...,m Il ) is an n — |m|-tuplc of in- 
tegers, each in a range of mt G {0, ±1, . . . , ±iV} and 
M(m) = X)"=i TOi - ^ main advantage of this method 
is that in contrast to H of Eq.©, H c g is constructed to 
conserve the total number of quanta Q at each order n. 
This is clear from the constraint M(m) = 0. Q conser- 
vation allows the SE of several quantities using the bare 
eigenstates of Ho like the ground state energy and the 
elementary triplet dispersion. This will be done in the 
next Sections. The amplitudes C(m) are rational num- 
bers computed from the flow equation method [lsl lifij . 
The C(m) table is available on |l7j . 
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GROUND STATE ENERGY AND ELEMENTARY 
TRIPLET DISPERSION 



FIG. 2: Ground state energy e g (a, b) — e g (a, 6 = 0) per tetra- 
hedron in the dimer phase. 



Now we discuss results for the ground state energy E g 
and the triplet dispersion ui(k) in the dimer phase as 
obtained from SE with respect to the inter-tetrahcdral 
coupling b. Q-conscrvation leads to 



E„ 



1-ffofflO), 



(5) 



where |0) = • • • <Si,£/2) and iJ ff is the effective 

Hamiltonian given by Eq.Q. Wrap-around of graphs up 
to length n will not occur if this matrix clement is evalu- 
ated on chains with n+1 tetrahcdral clusters and periodic 
boundary conditions (PBC). This is required by linkcd- 
cluster theorem and leads to SE's valid to 0(n) in 
the thermodynamic limit, i.e. for infinite-sized systems. 

We have evaluated E g up to 0(7). The ground state 
energy per tetrahedron e g in the dimer phase reads 



e g {a,b) = 



2 _ b l il 
" ~ ¥ + 108 
27311519 b 7 
1749600 + 5038848000' 



„ a 2b 

-2H 

2 3 

53273 b 6 



67 b 5 
1620 



(6) 



The first two terms in Eq.© correspond to the non- 
interacting energy of the Si state as in TablclJ For a = 
Eq. © is also a SE of the ground state energy per dimer 
of the dimerized spin— 1 chain. To the best of our knowl- 
edge an analytic expression of this has not been pub- 
lished previously. FigEl depicts the ground state energy 
which is a monotonously decreasing function of the inter- 
tetrahedral coupling b. We find, that plots of e g (a,b) to 
0(6) and 0(7) are indistinguishable on the scale of Fig|2 
which provides an estimate of convergence. Using e g {a, b) 
we will discuss the first order critical line b c (a c ) for the 
dimer-to-singlet transition in the next Section. Trivially, 
the non- interacting critical point is b c (l) = (table [IJ. 



To calculate the triplet dispersion we have to con- 
sider the subspace of exactly one 7[ j— type of excitation: 
\tj) = |6>i,i • ■ ■ Ti j . . . i.e. one-triplet 7ij on a 

tetrahedron j. By Q conservation the number of triplets 
is conserved, implying that H c r can only translate the 
triplets \tj) 



H eS \j) = ^ Cl \j 



(7) 



Due to translational invariance of our model, spatial 
Fourier transformation = y^E/LJZ • exp(ikj)\j), di- 
agonalizes Eq.Q leading to the triplet dispersion 



,{k) = (k\H cS \k)-Eg 

OO 

= c - E g + 2c j cos(fcj) . 

3=1 



(8) 



We have evaluated the coefficients Cj up to 0(7). From 
them, Eq.(Sj in the thermodynamic limit reads 
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(9) 

As for e g and since w(fc, 6) is independent of a, this result 
may also be interpreted as a SE for the triplet dispersion 
of a dimcrizcd spin— 1 chain. 

In FigO we compare our SE results (solid lines) for 
u>Qc,b) with other findings at an intermediate value of b 
which will be shown later to reside still within the dimer 
phase. The first excited eigenvalues obtained from exact 
diagonalization (ED) on a finite dimerized spin— 1 chain 
[l2| are displayed with dots. Results from bond-boson 
mean field theory (MFT) Il2l and linearized Holstein- 
Primakoff (LHF) methods TH |H Hi] are shown with 
gray and dashed lines, respectively. As we can observe 
in this Fig. the agreement between the SE's results and 
ED is very good, and the SE clearly improves on earlier 
uncontrolled approximations. Eq.© shows that there is 
a tendency to the closure of the gap: A = oj(k = 0, b) as 
b increases, which signals the transition to the Haldanc 
phase. This will be analyzed in detail in the next Section. 
Finally we note, that in the ID case our results extend 
to 0(7) an earlier SE which was obtained to 0(4) in ref. 
[lij . however for a 3D tetrahedral spin model. 



QUANTUM PHASE TRANSITIONS 

In this Section we consider the quantum critical points 
of our model. Fig0| summarizes the phase diagram of 
the tetrahedral cluster spin chain as determined from the 
present SE and other methods [lil fl3l ] . It displays the 
singlet product, dimcr and Haldane phases. These re- 
gions are separated by a first order critical line between 
the dimerized spin— 1 chain sector and the singlet prod- 
uct phase, as well as the second order critical line between 
the dimer and the Haldanc phases. 

We start by studying the dimer-to-singlet first order 
transition. The solid line represents our evaluation of 
the critical line b c (a c ) using SE. It has been obtained by 
numerical determination of 6(a) from e g (a, 6) — e g $ 2 = 0, 
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FIG. 3: Elementary triplet dispersion in the spin— 1 dimer 
phase for b = 0.3. Solid line: SE from Eq.@. Gray line, 
dashed line, and solid dots: bond-boson Mean Field Theory 
(MFT), bond-boson Linear Holstein Primakoff (LHP) approx- 



imation, and exact diagonalization (ED) from ref. 



3 Jar 



where e g s 2 = — 3a/2 is the exact ground state energy of 
the singlet-product phase. Obviously, there is an excel- 
lent agreement with ED results |12l l2"| (dotted line in 
FigEJ. Surprisingly this is true even for b- values beyond 
the dimer-Haldane transition, where the SE is not ex- 
pected to be valid anymore. From this we conclude that 
the ground state energy of the Haldane phase has lit- 
tle difference to that of an adiabatically continued dimcr 
phase. While the accuracy of approximate approaches, 
as eg. bond-operator MFT is hard to asses, it is 
interesting to see that the SE agrees well also with the 
latter approach at least up to the dimer-Haldane transi- 
tion (dashed line in Fig. 

The dimer-Haldane transition of the dimerized spin— 1 
chain has been analyzed by several authors. From the an- 
alytical point of view, Haldane and Affleck [26l Ul\ have 
performed a mapping onto the 0(3) nonlinear cr— model 
(NLSM) in the large S limit. They have shown that the 
topological angle 9 is given by 9 = 2ttS(1 — S) and that 
the system is gapless and described by a conformal field 
theory with SU(2) symmetry when 9/(2ir) is half odd 
integer. 8 is related to b by b = (1 — <5)/(l + S). There- 
fore, for S = 1 a gapless point was predicted to occur 
at b c = 1/3, i.e. 8 C = 1/2. While the prediction of a 
gapless point from the NLSM is expected to be correct, 
quantitative agreement for the critical point is not to be 
expected since 1 / 5*-correction may play a role and have 
not been analyzed to our knowledge. 

From the numerical point of view, the dimer-Haldane 
transition has been studied employing Density Ma- 
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FIG. 4: Quantum phase diagram of the tetrahedral cluster 
chain. Solid: first order dimer-to-singlet critical line b c (a c ) 
obtained from SE. Dashed and dots: b c (a c ) from bond-boson 
mean field theory (MFT) and exact diagonalization (ED) in 
the dimerized spin— 1 sector for N — 16 sites and PBC (ref. 
H3)- Horizontal solid line at b c £ [0.612, 0.614]: second order 
dimer-to-Haldane transition obtained by Dlog-Pade analysis 
of the triplet-gap closure of Eq.@. (See also Fig|SJ. 



trix Renormalization Grou p ( DMRG) j2j| |2jj, Quan- 
tum Monte Carlo (QMC) El, and ED techniques 
[l2l IsH IssT ] . These methods agree on a critical point at 
b c (S c ) ~ 0.6(0.25). SE has been applied to the dimer- 
ized spin— 1 chain [34| ■ From the SE of the ground state 
equal time structure factor and its second moment they 
obtained b c £ [0.56,0.64]. While the latter type of SE 
was based on ground state properties only, a SE for the 
elementary excitations and the triplet gap has not been 
given to our knowledge. 

Fig[S] shows the integrated (4, 3)-Dlog-Pade approxi- 
mant for the triplet gap from our SE at 0(7). As from 
Eq.©, the gap is located at k = 0, n. The axes have 
been scaled to allow for a comparison with two other 
approaches, i.e. ED on N = 16 sites with PBC in the 
dimerized spin— 1 chain sector an d DMRG results 
from ref. [28|. The latter two cover the complete range 
of b(S) G [0,1] ([1,0]). For b{5) < 0.4(> 0.4), there is 
a very good agreement among all three approaches dis- 
played. For greater(smallcr) values of b(S) finite size 
effects become evident in the ED data. The results 
from the SE however remains very close to the DMRG 
over the complete range of inter-tetrahedral couplings, 
up until the critical point. From the inset of Fig(5]one 
can observe that the SE slightly over(under)estimate the 
critical value of b c (S c ) as compared to the findings of 
DMRG. In fact, performing standard Dlog-Pade error 
analysis, by evaluating the scatter of the critical point 
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FIG. 5: Solid line: integrated Dlog-Pade (4, 3) approximant 
to the 0(7) SE, Eq.©, for the one-triplet gap A = uj(k = 
0,6). Axes have been scaled to allow for a comparison with 
other methods. Dots: ED on N — 16 sites with PBC in the 
dimerized spin— 1 sector [Q. Crosses: DMRG results from 
[28ll . Inset: zoom of the transition region. 



b c (8 c ) as obtained from different approximants, we get 
b c {S c ) 6 [0.612, 0.614] (e [0.239, 0.240]), whereas DMRG 
estimates b c (S c ) e [0.59, 0.61] (e [0.24, 0.26]) H|. The 
small difference may be due to non analytic corrections 

to a plain power-law behavior of A. Indeed, the gap 

2 1 

has been claimed to close as A ~ (A) 3 / logAlj , with 
A = 1 6— b c \/b c , in the vicinity of the critical point [13,I13- 
In case of such logarithmic corrections, deviations as 
those shown in the inset of Fig[5]arc likely to occur and 
even higher order SE would be required to improve the 
agreement very close the critical point. Further evidence 
for the relevance of logarithmic corrections stems from 
the critical exponent of v <E [0.98, 0.99] which we extract 
from the Dlog-Pade approximant which, as can be seen 
also from Fig. [5] varies almost linearly with b(8) close to 
the critical point. This is at variance with extrapolation 
of ED data, which predict a critical effective exponent 
v E [0.7, 0.8] (A ~ A") [13], and with the onset of ad- 
ditional curvature which can bc observed in the DMRG 
results close to b c (8 c ) in the inset of Fig|S] This issue 
should be addressed in future studies. 



CONCLUSIONS 

To summarize, we have studied zero temperature prop- 
erties of a tetrahedral cluster spin chain using a series 
expansion technique based on Wegner's flow equation 
method. Starting from the limit of isolated tetrahedra, 
we have obtained expansions up to 0(7) in the inter- 
tctrahcdral couplings for the ground state energy and 
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the dispersion of the elementary one-triplet excitations 
within the dimer phase. The ground state energy has 
been used to determine a first order quantum critical line 
which separates a dimerized spin— 1 chain sector from 
a singlet product phase. Our findings are in excellent 
agreement with those of ED. A second order critical line 
for a dimer-to-Haldane phase transition was obtained by 
analyzing the closure of the triplet gap. Again, very good 
agreement was found with the previous results from ED, 
DMRG and QMC. 
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